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We study the divergence of large-order perturbation theory in the worldline expression 
for the two-loop Euler-Heisenberg QED effective Lagrangian in a constant magnetic field. 
The leading rate of divergence is identical, up to an overall factor, to that of the one-loop 
case. From this we deduce, using Borel summation techniques, that the leading behaviour 
of the imaginary part of the two-loop effective Lagrangian for a constant E field, giving the 
pair-production rate, is proportional to the one-loop result. This also serves as a test of the 
mass renormalization, and confirms the earlier analysis by Ritus. 
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I. INTRODUCTION 

Euler and Heisenberg and many others since computed the exact renormalized one-loop 

I ' QED effective action for electrons in a uniform electromagnetic field background. When the background is 
. purely that of a static magnetic field, the effective action is minus the effective energy of the electrons in 
that background. When the background is purely that of a uniform electric field, the effective action has an 
imaginary part which determines the pair- production rate of electron-positron pairs from vacuum [Qj2|^,^,||,|| . 
In this paper we consider the two-loop Euler-Heisenberg effective action, and we show how the divergence 
of the perturbative expression for the effective action with a uniform magnetic background is related to the 
non-perturbative imaginary part of the effective action with a uniform electric background. The two-loop 
Q\ , Euler-Heisenberg effective Lagrangian, describing the effect of a single photon exchange in the electron loop, 
was first calculated by Ritus 0, and later recalculated by Dittrich and Reuter H for the magnetic field 
case. In both cases the proper-time method and the exact Dirac propagator in the uniform field were used. 
More recently the magnetic field computation was repeated in || using the more convenient 'worldline' 
formalism 0,0. This calculation revealed that the previous results by Ritus and Dittrich/Reuter were 
actually incompatible, and differed precisely by a finite electron mass renormalization. This prompted yet 
another recalculation of this quantity in the worldline formalism [ [l2| , now using dimensional regularisation 
instead of a proper-time cutoff as had been used in the previous calculations. That calculation confirmed 
the correctness of Ritus's result, and conversely showed that the final result given by §| was not expressed 
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in terms of the physical electron mass. As part of our analysis here, we show how this finite difference in 
the mass renormalization affects the large-order behaviour of perturbation theory, and how this affects the 
leading contribution to the imaginary part of the effective action in the electric field case. 

For a uniform magnetic background, of strength B, the one-loop effective Lagrangian has a simple "proper- 
time" integral representation : 

The i term is a subtraction of the zero field (B — 0) effective action, while the | subtraction corresponds to 
a logarithmically divergent charge renormalization. For a given strength B, this integral can be evaluated 
numerically || . Alternatively, we can make contact with a perturbative evaluation of the one-loop effective 
action by making an asymptotic expansion of the integral in the weak field limit - i.e., for small values of 
the dimensionless parameter 
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Here the 2?2n are Bernoulli numbers Each term in this expansion of £W is associated with a one-fermion- 
loop Feynman diagram. Note that only even powers of eB appear, as expected due to charge conjugation 
invariance (Furry 's theorem). The divergent 0(e 2 ) self-energy term is not included as it contributes to the 
bare Lagrangian by charge renormalization. 



The expansion (|2j) is the prototypical "effective field theory" effective Lagrangian 14 1, where the low 



energy effective Lagrangian, for energies well below the fermion mass scale m, is expanded as 

Q( n ) 

£ = m iJ2a n — (3) 

n 

with C)( n ' being an operator of dimension n. For QED in a uniform background, the higher dimensional 
operators are formed from powers of Lorentz invariant combinations of the uniform field strength F^ u . 
For a uniform magnetic background this simply means even powers of B 7 as in (|^). Note that the 'low 
energy' condition here means that the cyclotron energy — is well below the fermion mass scale m; in other 
words, ^ < 1. 

The Euler-Heisenberg Lagrangian encodes the information on the low-energy limit of the one-loop N - 
photon amplitudes in a way which is highly convenient for the derivation of various nonlinear QED effects 
such as vacuum birefringence (see, e. g., [[l5| and refs. therein) or photon splittin g ]l^,[l7t . The experimental 
observation of vacuum birefringence is presently attempted by laser experiments Jl8| , [l9| . There is also recent 
experimental evidence for vacuum effects in pair production with strong laser electric fields p0| . 

The one-loop Euler-Heisenberg perturbative effective action (||) is not a convergent series. The one- loop 
expansion coefficients in (||) alternate in sign [since sign(£>2„) = (— l) n+1 ]) but grow factorially in magnitude 
(see also Table 1): 

(1) _ 2 2 ™£>2n+4 

fl " ~ ~(2?i + 4)(2n + 3)(2n + 2) 

1 T(2n + 2) /ll 



~ ( n 2n y 1 + 2 2«+4 + 3 2«+4 +"•) ( 4 ) 

So the perturbative expansion (^) is a divergent series. This divergent behaviour is not a bad thing; it is 
completely analogous to generic behaviour that is well known in perturbation theory in both quantum field 
theory and quantum mechanics pl| . For example, Dyson ]22| argued physically that QED perturbation 
theory is non-analytic, and therefore presumably divergent, as an expansion in the fine structure constant 
a, because the theory is unstable when a is negative. As is well known, the divergence of high orders of 
perturbation theory can be used to extract information about non-perturbative decay and tunneling rates, 
thereby providing a bridge between perturbative and non-perturbative physics pl[] . It has been argued p3| , 
based on the behaviour of the one- loop Euler-Heisenberg effective Lagrangian (H), that the effective field 
theory expansion (|^) is generically divergent. Here we consider this question at the two- loop level. 

We stress that for energies well below the scale set by the fermion mass m, the divergent nature of the 
effective Lagrangian is not important, as the first few terms in the series (||) provide an accurate approxima- 
tion. However, the divergence properties do become important when the external energy scale approaches 
the fermion mass scale m. The divergence is also the key to understanding how non-perturbative imaginary 
contributions to the effective action arise from real perturbation theory. 



II. BOREL ANALYSIS OF THE ONE-LOOP EULER-HEISENBERG EFFECTIVE LAGRANGIAN 



To begin, we review very briefly some basics of Borel summation |2j|,|2jj|26|,|27| . Consider an asymptotic 
series expansion of some function f(g) 
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where <? — > + is a small dimensionless perturbation expansion parameter. In an extremely broad range of 
physics applications |2l| one finds that perturbation theory leads not to a convergent series but to a divergent 
series in which the expansion coefficients a n have leading large-order behaviour 



a„~ (-l)VT(/xn + i/) 
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for some real constants p, p, > 0, and z/. When p > 0, the perturbative expansion coefficients a n alternate 
in sign and their magnitude grows factorially, just as in the Euler-Heisenberg case (f|). Borel summation is 
a useful approach to this case of a divergent, but alternating series. Non-alternating series must be treated 
somewhat differently. 

To motivate the Borel approach, consider the classic example : a n = (—l) n p n nl, and p > 0. The series 
(H) is clearly divergent for any value of the expansion parameter g. Write 
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where we have formally interchanged the order of summation and integration. The final integral, which is 
convergent for all g > 0, is defined to be the sum of the divergent series. To be more precise , the formula 

(0) should be read backwards: for g — > + , we can use Laplace's method to make an asymptotic expansion 
of the integral, and we obtain the asymptotic series in (^) with expansion coefficients a n — (— \) n p n n\. 

For a non-alternating series, such as a n = p n n\, we need f(—g). The Borel integral (Q) is [p4|]25| l an 
analytic function of g in the cut g plane: |arg(g)| < n. So a dispersion relation (using the discontinuity 
across the cut along the negative g axis) can be used to define the imaginary part of f(g) for negative values 
of the expansion parameter: 
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The imaginary contribution (||) is non-perturbative (it clearly does not have an expansion in positive powers 
of g) and has important physical consequences. Note that (^) is consistent with a principal parts prescription 
for the pole that appears on the s > axis if we make the formal manipulations as in (0): 
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Similar formal arguments can be applied to the case when the expansion coefficients have leading behaviour 
J). Then the leading Borel approximation is 
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For the corresponding non-alternating case, when g is negative, the leading imaginary contribution is 
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Note the separate meanings of the parameters p, p, and v that appear in the formula (|^) for the leading large- 
order growth of the expansion coefficients. The constant p clearly combines with g as an effective expansion 
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parameter. The power of the exponent in ( |11| ) is determined by p, while the power of the prefactor in ( |11| ) 
is determined by the ratio ^. 

It must be stressed that these formulas (|l^) and ([n]) are formal, being based on assumed analyticity 
properties of the function f(g). The Borel dispersion relations could be complicated by the appearance of 
additional poles and/or cuts in the complex g plane, signalling new physics p7| . In certain special cases these 
analyticity assumptions can be tested rigorously, but we have in mind the situation in which one is confronted 
with the expansion coefficients a n of a perturbative expansion, without corresponding information about the 
function that this series is supposed to represent. This is a common circumstance in physical applications of 
perturbation theory. For example, Borel techniques have recently been used to study the divergence of the 
derivative expansion for QED effective actions in inhomogeneous backgrounds ]2§| ] . 

Returning to the Euler-Heisenberg effective Lagrangian, the question of whether the perturbative expan- 
sion is alternating or non-alternating is directly relevant. For a uniform magnetic background, the one-loop 
Euler-Heisenberg series (||) is precisely of the form (|^) with g = (f^) 2 - Moreover, from (^) the expansion 

coefficients have leading large-order behaviour of the form (^), with p = and p = v = 2. In fact, 
taking into account the sub-leading corrections indicated in (I), the proper-time integral representation (Jlj) 
is precisely the Borel sum, using (|l0|), of the divergent series (2) [ p8[ . For a uniform electric background, the 
only difference perturbatively is that B 2 is replaced by — E 2 ; that is, g = {^) 2 is replaced by — g = —(f^) 2 - 
So the perturbative one-loop Euler-Heisenberg series (Q) becomes non-alternating. Then from (^), with 
p = and p = v = 2, we immediately deduce the leading behaviour of the imaginary part of the one-loop 
Euler-Heisenberg effective Lagrangian: 
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This imaginary part has direct physical significance - it gives half the electron-positron pair production rate 
in the uniform electric field E ||||J{|. Actually, since we also know the sub- leading corrections (|]) to the 



leading large-order behaviour of the expansion coefficients a n 
the leading behaviour in (|L2]): 
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we can apply (nil) successively to go beyond 
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This is Schwinger's classic result || for the imaginary part of the one-loop effective Lagrangian in a uniform 
electric field E. To elucidate the physical meaning of the individual terms of this series it is useful to employ 
the following alternative representation due to Nikishov B , 
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Here ftp is the mean number of pairs produced by the field in the state with given momentum p and spin 
projection r. An expansion of the logarithm in rip and term- by-term integration leads back to Schwinger's 
formula (|l3|). Thus the leading term in this formula can be interpreted as the mean number rip of pairs in 
the unit 4-volume VT, while the higher (k > 2) terms describe the coherent creation of k pairs. 

Pair creation can occur for any value of the electric field strength, though due to the exponential suppression 
factors one is presently still far away from being able to observe spontaneous pair creation by macroscopic 
fields in the laboratory. However, it can be arranged for electrons traversing the focus of a terawatt laser 
to see a critical field in their rest frame. This has recently led to the first observation of pair creation in a 
process involving only real photons pp[ |. 

For the one-loop Euler-Heisenberg QED effective Lagrangian, this large-order perturbation theory analysis 
is greatly simplified by the fact that we know the exact formula (|J) for the expansion coefficients al . This 
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will not be the case below, when we discuss the two-loop Euler-Heisenberg effective Lagrangian. So, for the 
sake of numerical comparison, we compare the exact one-loop coefficients affl with their leading large-order 
behaviour. The coefficients are listed in Table 1 up to n — 15. Since the growth is fast, it is convenient to 
compare the logarithms, as is done in Figure 1. With 16 terms it is straightforward to fit the the values 
of p, /i and v appearing in (^|); moreover, there is sufficient accuracy to fit the overall coefficient g^pj. In 

Figure 1 we plot = log |a!^|, and = log[r(2n + 2)/(87r 2n+4 )]. The agreement is spectacular, already 
for n = 0. Indeed, on this scale the two plots are indistinguishable. To go beyond the leading large-order 
behaviour, we plot the difference - C$P . This can be fitted to the correct form log(l + ga^ra ) ~ 2^"+ 4 ' 
with remarkable accuracy for n > 2, as illustrated in Figure 2. 

TABLE I. The one-loop Euler-Heisenberg coefficients afi from (Q) and their magnitudes |o^|. The last two 
columns list the calculated two-loop Euler-Heisenberg coefficients a! 2 ' in (pi]), and their magnitudes Note that 

both and a£ 2 ^ alternate in sign and grow factorially in magnitude. 
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FIG. 1. Plots of An' = log(|ak ; |) and C n ' = log[2r(2n + 2)/(87r 2n+4 )] up to n = 15. The two plots are 
indistinguishable on this scale. 




FIG. 2. Plot (solid line) of the difference A n 1} - C ( n L> between A 



log(|on"^|) and the leading large order fit 



Ci 1) = log[r(2n + 2)/(87r 2n+4 )] from (g), up to n = 15. The dashed line shows the fit of this difference to the function 
log(l + 2 il+i ), which is excellent already for n > 2. 



III. THE TWO-LOOP EULER-HEISENBERG EFFECTIVE LAGRANGIAN 

We now turn to the two-loop Euler-Heisenberg effective Lagrangian, describing the effect of a single 
photon exchange in the electron loop. This quantity was first studied by Ritus 0. Using the exact electron 
propagator in a constant field found by Fock |29| and Schwinger ||, and a proper-time cutoff as the UV 
regulator, he obtained the on-shell renormalized in terms of a certain two-parameter integral. From this 
integral the imaginary part of the Lagrangian was then extracted by a painstaking analysis of its analyticity 
properties, yielding a representation analogous to Schwinger's one-loop formula ( |l3| ) pof . Adding up the 
one-loop and the two-loop contributions to the imaginary part the result reads 
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where (3 



For the function (/3) the following small (3 - expansion was obtained in p0| 
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According to M the physical interpretation of the individual terms in the series in terms of coherent multipair 
creation can be carried over to the two-loop level. This requires one to absorb the term linear in a into the 
exponential factor, rewriting 
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The individual terms in the expansion ( |15|) of Kk{P) arc then to be absorbed into the mass shift m* — m. 
For the lowest order terms in this expansion, those given in ([l5|), the physical meaning of the corresponding 
mass shifts in terms of the coherent pair production picture is discussed in [Q. For example, the leading 
"1" in the expansion of Ki([3) after exponentiation yields a mass shift that can be identified as the classical 
change in mass caused for one particle in a created pair by the acceleration due to its partner. 

Assuming this exponentiation to work one can, of course, obtain some partial information on the higher- 
loop corrections to the imaginary part. More remarkably, since the above physical interpretation requires the 
mass appearing in the exponent to be the physical one, it allows one to determine the physical renormalized 
mass from an inspection of the renormalized Lagrange function alone, rather than by a calculation of the 
(lower order) electron self energy. 

Following the pioneering work by Ritus and his collaborators, a first recalculation of the Euler-Heisenberg 
Lagrangian was done by Dittrich and Reuter || for the magnetic field case. The more recent recalculation 
in |^] showed that the two previous results were actually incompatible, and differed precisely by a finite 
electron mass renormalisation. All three calculations had been done using a proper-time cutoff rather than 
dimensional regularisation. This cutoff leads to relatively simple integrals, but due to its non-universality 
makes it, at the two-loop level, already non-trivial to determine the physical renormalized electron mass. A 
fourth calculation of this quantity ]l2[ |, now using dimensional regularisation, yielded complete agreement 
with Ritus's result after a perturbative expansion of both results in powers of the B field had been performed. 
In the following we will push the same calculation to 0(B 34 ), and analyze the rate of growth of the expansion 
coefficients. 



IV. BOREL ANALYSIS OF THE TWO-LOOP EULER-HEISENBERG EFFECTIVE 

LAGRANGIAN 



The world-line expression for the two-loop on-shell renormalized Euler-Heisenberg effective Lagrangian is 
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Here a w is the fine-structure constant, and 7 = 0.5772... is Euler's constant. The function L(s,u) 
appearing in the integrand of the first term in ( jl8| ) is defined by the following relations: 
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(19) 
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The second term in the two-loop expression (|18j) is generated by the one-loop electron mass renormalisation, 
which at the two-loop level becomes necessary in addition to the photon wave function renormalisation. For 
this mass renormalization term we have found the following exact expansion: 
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Here ii?2n are the Bernoulli numbers, and ip(x) = jr^j is the digamma function [Q. We have not succeeded 
in finding a closed-form expression for the expansion coefficients arising from the expansion of the double 
integral in (|l8|), although we suspect that one may exist. Instead, we have made an algebraic expansion of 
this integral, using MATHEMATICA and MAPLE. When combined with the exact expansion (|o|) of the 
mass renormalization term we obtain an expansion of the form: 
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The expansion coefficients a n are listed in Table 1, up to n = 15. Note that those coefficients are in 
some sense less universal than their one-loop counterparts, since they depend on the one-loop normalization 
condition imposed on the renormalized electron mass. 

Several comments are in order. First, the two-loop expansion coefficients a n alternate in sign, just as 
in the one- loop magnetic background case (^). Second, the magnitude \a^\ is clearly growing factorially 
fast with n. Thus, the two-loop Euler-Heisenberg series ( |2l] ) is a divergent series, as is the one-loop Euler- 
Heisenberg series (||) . Note also that for each series the smallest magnitude coefficient is reached already for 
n=l, after which the coefficients begin to increase rapidly in magnitude. 

To extract the leading large-n growth of la^l we fit to the form in (|^). Once again, it is convenient 
to work with the logarithm D n — log(|a„|) since the growth is so rapid. It is relatively straightforward to 
find that fi = v = 2 and p = -\. It is more difficult to fit the overall coefficient, but if we assume this is a 
simple power of ir then our best fit for the leading large-order growth of the two-loop expansion coefficients 
in (pit) is: 
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This leading fit is displayed in Figure 3, in terms of An = log(|a„ |) . The fit is not as good as the one-loop 
fit shown in Figure 1, but it is still very good. 

Note the remarkable similarity of the leading large-order growth ( |22[ ) of the two-loop expansion coefficients 
to the leading large-order growth of the one-loop expansion coefficients in (Q). The only difference is the 
overall coefficient. The parameters p, p and v in the general form (^) are identical. Using the Borel technique 
to relate this leading growth rate to the leading non-perturbative imaginary part of the effective Lagrangian 
in a uniform electric field E, we deduce that the two-loop leading imaginary part is proportional to the 
one-loop leading imaginary part (|l2"|). In fact, from ( jll| ) and (22), we find the leading contribution 
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FIG. 3. Plots of A ( n } = log(|ai 2) |) (solid line), and C*i 2) = log[16r(2n + 2)/7r 2n+2 ] (dashed line), up to n = 15. 
While the agreement is not as good as for the one-loop case shown in Figure 1, the fit is still very good, even at low 
orders in n. 



This has exactly the same dependence on the electric field E as the one-loop case. So to two-loop order 
the leading non-perturbative behaviour of the imaginary part of the effective Lagrangian is: 



(24) 



Im (£«+£t 2 
This agrees with the leading term of Ritus's formula ( [To] ) 
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FIG. 4. Plot (solid line) of the difference — C)?' between An' — log(|a^ z; |) and the leading large order fit from 
@, C ( n ] = log[16r(2rz + 2)/7r 2n+2 ], up to n = 15. The dashed line represents the fit in (||). 

To go beyond this leading term we need to look at corrections to the leading behaviour in (^3|) . In Figure 
4 we plot the difference of the logarithms, and we see that the n dependence is much more gentle than the 
rapid fall-off found in the one-loop case, which was plotted in Figure 2. In fact, from the terms up to n = 15, 
we obtain the fit 
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This is a considerably weaker n dependence than is found for the first correction in the one- loop case (Q). 
This means that in the two-loop case the dominant corrections are to the prefactor in the leading behaviour 
(p3|). This is in contrast to the one- loop case ( [l3| ) where the first correction to the leading behaviour is 
exponentially suppressed. Indeed, applying the Borel relations, the correction term (Gq) leads to 
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We emphasize that the fit in ( j25| ) is based on a simple fit to the first 16 two-loop coefficients. Nevertheless, 
the structure of ( |26] ) conforms already to the form of Ritus's expansion in eq. (|l5|) . It would be interesting 
to probe this correction term in more detail by a further study of the analyticity properties of the integral 
representations [|]],||,|l2| of the two-loop Euler-Heisenberg effective Lagrangian, or by looking at still higher 
orders in perturbation theory 



V. CONCLUDING REMARKS 



Our analysis also permits us to study the dependence of (|23|), the leading non-perturbative imaginary 
contribution to the effective Lagrangian, on the electron mass renormalisation. In the world-line expression 
(|l8|) for the two-loop Euler-Heisenberg effective Lagrangian, a finite change of the renormalised electron mass 



would amount to an arbitrary change of the constant 5 — 7 appearing in the second bracket of the second 
term. For example, in |9| it had been shown that the renormalised Lagrangian obtained by [§] differs from 
( |l8| ) precisely by a replacement of | by |. A separate study of the contributions of the first and second 
term in ( |l8| ) shows that, due to cancellations between both terms, the leading large-n growth of their sum 
is smaller than for each term separately. However this property holds true only if the renormalised electron 
mass is the physical one. 




FIG. 5. Plot (dashed line) of the correct leading behaviour, C { „ 2) = log[16r(2n + 2)/7r 2n+2 ], of the logarithm of 
the magnitude of the two- loop expansion coefficients; as compared to (solid line) = log(jfi^|), obtained from 
the two- loop coefficients off extracted from the result of One can clearly see the divergence between these two 
curves at large n, indicating a difference in the leading growth of the coefficients. 

For definiteness, in Figure 5 we compare the correct leading growth ( p2|) of the two- loop coefficients with 
the coefficients obtained by expanding out the two- loop result of [||. The difference in the leading large-n 
growths is obvious. Thus, in agreement with the analysis of Q, we find that if and only if expressed in terms 
of the true electron mass will the imaginary part of the renormalised two-loop Lagrangian show the same 
exponential suppression factor exp[— irm 2 /(eE)] as for the one-loop Lagrangian. 

To summarize, we have constructed the imaginary part of the two-loop QED Euler-Heisenberg Lagrangian 
in a constant electric field by a computer based calculation of its weak field expansion together with Borel 
summation techniques. The knowledge of the first 16 coefficients has turned out to be sufficient to verify 



the structure of the leading (k = 1) term in Ritus's eq.(15), and to obtain a numerical value for the first 
0{^JeE /m 2 ) correction contained in that formula. The method used here is significantly simpler than the 
one in 0,^O), where the imaginary part was obtained by an analysis of the analyticity properties of the 
two-loop parameter integrals. 

In particular, we have seen that the large order behaviour of the two-loop coefficients is the same (up to an 
overall constant factor) as the one-loop case. This means that the leading contribution to the imaginary part 
of the two-loop effective Lagrangian has the same form as in the one-loop case. This gives a new perspective 
to Ritus's arguments j7j that the true renormalized electron mass m is such that the leading exponential 
factor in the pair production rate is exp[— irm 2 / (eE)]. Since those arguments pertain to arbitrary loop orders, 
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and the leading exponential factor is directly related to the leading growth rate in the weak field expansion, 
they also lead one to expect that the Euler-Heisenberg Lagrangian may be amenable to this type of Borel 
analysis at any fixed loop order in perturbation theory. 



VI. ACKNOWLEDGEMENTS: 



C. S. would like to thank R. Stora for discussions. This work has been supported in part (GD) by the 
U.S. Department of Energy grant DE-FG02-92ER40716.00. 



[1] W. Heisenberg and H. Euler, "Folgerungen aus der Diracschen Theorie des Positrons", Z. Phys. 98, 714 (1936). 
[2] V. Weisskopf, "Uber die Elektrodynamik des Vakuums auf Grund der Quantentheorie des Elektrons", Kong. Dans. 

Vid. Selsk. Math-fys. Medd. XIV No. 6 (1936), reprinted in Quantum Electrodynamics, J. Schwinger (Ed.) (Dover, 

New York, 1958). 

[3] J. Schwinger, "On Gauge Invariance and Vacuum Polarization", Phys. Rev. 82, 664 (1951). 
[4] A. Nikishov, "Pair production by a Constant External Field", Sov. Phys. JETP 30, 660 (1970). 
[5] W. Dittrich, "One loop effective potentials in QED", J. Phys. A9, 1171 (1976). 
[6] W. Greiner and J. Reinhardt, Quantum, Electrodynamics (Springer, 1992). 

[7] V. I. Ritus, "Lagrangian of an intense electromagnetic field and quantum electrodynamics at short distances", Sov. 
Phys. JETP 42, 774 (1976); V. I. Ritus, "Connection between strong-field quantum elctrodynamics with short- 
distance quantum electrodyanmics" , Sov. Phys. JETP 46, 423 (1978); V. I. Ritus, "Effective Lagrange Function 
of Intense Electromagnetic Field in QED", Proc. of Frontier Tests o f Quantum Electrodynamics and Physics of 
the Vacuum, Sandansky 1998, 1998 Heron Press, 11, |hep-th/98 12124 . 

W. Dittrich and M. Reuter, Effective Lagrangians in Quantum Electrodynamics, (Springer, 1985). 
M. Reuter, M. G Schmidt and C. Schubert, "Constant Exte rnal Fields in Gauge Theory and the spin 0, 1/2, 1 
Path Integrals", Ann. Phys. 259, 313 (1997), frep-th/9610191 . 

M. J. Strassler , "Field theory without Feynman diagrams: one loop effective actions", Nucl. Phys. B 385, 145 
(1992). 

C. Schubert, "An Introduction to the Worldline Technique for Quantum Field Theor y Calculations", Lectures at 



XXXVI Cracow School of Theoretical Physics, Acta Phys. Polon. B27, 3965 (1996), bep-th/9610108 



D. Fliegner, M. Reuter, MG. Schmidt and C. Schubert, " The two loop Eul er-Heisenberg Lagrangian in dimensional 
renormalization" , Theor. Math. Phys. 113, 1442 (1997), jhep-th/9704194 . 

I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, Academic Press, New York, 1979. 
J. Donoghue, E. Golowich and B. Holstein, Dynamics of the Standard Model (Cambridge Univ. Press, 1992). 
W. Dittrich and H. Gies, "Vacuum Birefringence in Strong Magnetic Fields", Workshop Sandansky 1998, loc. cit., 
29, |hep-ph/9806417 . 

Z. Bialynicka-Birula and I. Bialynicka-Birula, Phys. Rev. D 2 (1970) 2341. 

SL. Adler, "Photon Splitting and Photon Dispersion in a Strong Magnetic Field", Ann. Phys. NY 67 (1971) 599. 
R. Pengo et al., "Magnetic Birefringence of Vacuum: the PVLAS Experiment", Workshop Sandansky 1998, loc. 
cit., 59. 

S.A. Lee et al., "Measurement of the Magnetically-induced QED Birefringence of the Vacuum and an Improved 
Laboratory Search for Light Pseudoscalars" , FERMILAB-PROPOSAL-P-877A, 1998. 

D.L. Burke et al., "Positron Production in Multiphoton Light-by-Light Scattering", Phys. Rev. Lett. 79 (1997) 
1626. 

For an extensive review, see: J. C. Le Guillou and J. Zinn- Justin (Eds.), Large-Order Behaviour of Perturbation 
Theory, (North Holland, Amsterdam, 1990). 

F. J. Dyson, "Divergence of Perturbation Theory in Quantum Electrodynamics", Phys. Rev. 85, 631 (1952). 
A. Zhitnitsky, "Is an Effective Lagrangian a Convergent Series?", Phys. Rev. D 54, 5148 (1996). 

G. Hardy, Divergent Series (Oxford Univ. Press, 1949). 



11 



[25] 

[26] 
[27] 



[28] 



C. M. Bender and S. A. Orszag, Advanced Mathematical Methods for Scientists and Engineers (McGraw-Hill, New 
York, 1978). 

J. Zinn- Justin, Quantum Field Theory and Critical Phenomena (Oxford Univ. Press, 1997). 

See for example: G. 't Hooft, in The whys of subnuclear physics, Erice Lectures 1977, A. Zichichi (Ed.), (Plenum 
Press, 1979); V . I. Zakharov, "Renormalons as a Bridge bet ween Perturbative and Nonperturbative Physics", 
hep-ph/9802416| M. Beneke, "Renormalons", |hep-ph/9807443| . 



G. Dunne and T. Hall, "Borel Summation of the Derivative Expansion and Effective Actions", hep-th/9902064 
Phys. Rev. D in press. 
V. Fock, Phys. Z. Sowjet. 12 (1937) 404. 

S.L. Lebedev and V.I. Ritus, "Virial representation of the imaginary part of the Lagrange function of the electro- 
magnetic field", Zh. Eksp. Teor. Fiz. 86 (1984) 408, [Sov. Phys. JETP 59 (1984) 237]. 



12 



